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ABSTRACT 



We review some universal features of the colour flux tube of gauge theories in the 
confining phase predicted by the infrared conformal limit of the underlying string 
theory. In particular we discuss shape effects in Wilson loops and rederive in a 
general way the logarithmic growth of the mean square width of the flux tube as 
a function of the interquark separation. Recent data on 3-D Z2 gauge theory, 
combined with high precision data on the interface physics of the 3-D Ising model 
fit nicely to this behaviour over a range of more than two orders of magnitude. 



1. Introduction 

1.1. The String Picture of the Flux Tube 

Experience on lattice gauge theory has shown that the relevant degrees of freedom 
of whatever gauge system in the confined phase are concentrated, in the infrared 
region, inside string-like flux tubes connecting among them the sources of the gauge 
held. This suggests, according to an old conjecture 1 ' 2 , that the vacuum expectation 
values of large Wilson loops W(C) } once the integration on the irrelevant degrees of 
freedom has been performed, are given by an effective string theory describing the 
dynamics of the flux tube as follows 



where E is a surface bounded by the loop C and S is a (largely unknown) string 
action. Within SU(N) gauge theories there is an internal support of this idea: the 
1/N expansion of weak coupling perturbation theory can be interpreted, according 
to 't Hooft analysis 2 , as a topological expansion of a string theory in the number of 
handles of E in which 1/N acts as the string coupling g. For large N the vacuum 
expectation value (W(C)) is dominated by planar Feynman graphs. In the dual 
language of strings this implies that the sum in Eq. I is saturated by surfaces with 
no handles and the evaluation of (W(C)) is greatly simplified. 

It has been recently shown 3 that in two-dimensional QCD on a Riemann surface 
the coefficients of the 1/N expansion have a simple, exact interpretation in terms of 
a string theory. 
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1.2. Condensation of Handles 

For finite N, or for discrete gauge groups, the string enters a non- perturbative 
regime where the creation of handles becomes an important process. For this reason 
a stringy description of such gauge systems, which are the most interesting from the 
physical point of view, is often considered hopeless. 

There is however an old argument 4 , corroborated by recent numerical simula- 
tions 5 ' 6 , which suggests that this non- perturbative regime has actually a simple 
description. Indeed it has been shown long ago 4 that in critical string theory the 
sum of the divergent part of the handle insertions in any string amplitude can be 
reabsorbed in the renormalization of the string tension a and of the string coupling 
g. An indirect check of this property comes from the simplest gauge system, namely 
the three-dimensional Z 2 gauge model. Here, using the duality relation with the 
Ising model, it is possible to describe the properties of the effective string in terms of 
the fluid interfaces between domains of opposite magnetization. It is also possible to 
study the topological properties of these surfaces and it turns out 5 that the number of 
handles is very large and proportional to the area of the surface. Such a condensation 
of handles seems an ubiquitous phenomenon of the physics of random surfaces with 
fluctuating topology 6 ' 7 ; its meaning is rather simple: the mean handle is microscopic 
and its size I does not depend on that of the surface. As a consequence, at a scale 
larger than I there are very few visible handles, so the effective string coupling g } being 
proportional to the amplitude of handle formation, is vanishing, like the perturbative 
coupling in the large A" gauge systems. 

The above argument strongly suggests that at large distances the functional in- 
tegral given in Eq. 1 is dominated, like in the 1/N expansion, by surfaces with no 
handles, so the effective string description of the infrared region should greatly sim- 
plify. 

1.3. a' Expansion and the Universality Class of the Flux Tube 

Linear confinement implies that the action S is proportional, in the infrared region, 
to the area of the surface E swept by the flux tube, then the string action should be 
well approximated by the Nambu-Goto action, even if this simple string picture does 
not work at short distances 8 . 

A further simplification of the string action comes from the so called a! expansion, 
where a! = l/2ircr : for large rectangular Wilson loops of size L\ X L 2 one can expand 
the action S in the natural adimensional parameter l/(crLiL 2 ) } obtaining 

S = °UU + \ f Ll d^ f L2 d£ 2 EW(£i,£ 2 )) 2 + O {-^-A , (2) 
2 Jo Jo ~~ \uL\L 2 J 

where hi are the two-dimensional fields describing the transverse displacements of the 
flux tube. Thus the effective string theory reduces in the infrared limit to a massless 
gaussian model. 

The above conclusion can be actually obtained also by more general considera- 
tions. Indeed strong coupling expansions of (W(C)) in whatever gauge theory 



show that the flux tube undergoes a transition towards a rough phase for (3 larger 
than a given fi r (the roughening point). It is widely believed that this transition is 
described by the universality class of Kosterlitz and Thouless 9 (KT). Thus the ef- 
fective string action, even if it is substantially unknown, should belong to the same 
universality class of the XY model. The renormalization group equations of the KT 
universality class can be simply expressed in terms of two couplings: 

x(t) = -y(tf , (3) 
y(t) = -x(t)y(t) , (4) 

where t is the RG scale parameter, x = irT — 2, (T is the temperature of the XY 
model) and y is the vortex fugacity. The inverse temperature can be identified in 
the infrared limit with the surface stiffness k, which in the continuum limit coincides 
with the string tension. These RG equations describe the transition from a smooth 
phase (with free vortices) to a rough phase where the effective vortex fugacity flows 
to zero and then the the theory becomes at large scale a massless free field theory. 

So, we have reached a conclusion similar, but more general than that obtained by 
the a' expansion of the Nambu-Goto action: for large enough Wilson loops it is not 
necessary to know the specific form of the string action because it flows towards a 
free massless theory. This limit action is of course critical in the whole rough phase, 
then it produces universal (i.e. gauge group independent) finite size effects. Moreover, 
because this theory is also infrared divergent, such finite size effects are expected to be 
rather strong. We shall discuss two of these universal effects: the shape dependence 
of the Wilson loops and the growth of the mean width of the flux tube. 



2. Shape Effects in the Wilson Loops 

2.1. The Functional Form of the Quantum Fluctuations 

The contribution of the first few terms of the a! expansion to the vacuum expec- 
tation values of the Wilson loops can be calculated exactly 10 . Here we consider only 
the contribution of the gaussian term, which is the one associated to the universal 
effects. 

For large enough Wilson loops we can parametrize (W(L\, L 2 )} as follows 

- log(W(L 1 ,L 2 )) = aL x L 2 + p(Lr + L 2 ) + c + q{L u L 2 ) + O(a') (5) 

where the string tension a and the the two parameters c and p are the only free 
parameters of this equation; q(L\, L 2 ) is the universal contribution due to the gaussian 
functional integral and is a know function of the shape (L\, L 2 ) which does not depend 
on the gauge system, but only on the choice of the boundary conditions. For instance, 
for fixed boundary conditions we have 

q{ x i y) = lo s I ^ J , (6) 



where the Dededekind eta function is given by 
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The older fits of Wilson loop data put arbitrarily q(x } y) = with unacceptable 
X 2 values. Insertion of the quantum correction given in Eq. 6 improves considerably 
the fits, even if discrepancies have been observed for various gauge systems both in 
D = 3 and D = 4 dimensions n > 12 . 

It has been suggested that these discrepancies might be due to the choice of 
fixed b.c. 13 . Actually the gaussian model is a member of a one-parameter family of 
conformal field theories ( those with central charge c = 1) having the same local 
behaviour, but different topological modes (winding modes around a circle of radius 
r which is the parameter of this family). If the free theory describing the infrared 
limit of the effective string is allowed to have such topological modes, it is possible to 
introduce new kinds of boundary conditions which give a much better description of 
the Wilson loop. Roughly speaking, such a topological modification of the gaussian 
model allows one to treat in a different manner quarks line and antiquarks, while in 
the ordinary gaussian model there is no way to distinguish a flux tube connecting two 
quarks from that connecting a quark-antiquark pair, which from the point of view of 
the physics of the gauge system are of course completely different. 

2.2. Fluid Interfaces 

In order to test more accurately the above universal asymptotic behaviour it would 
be convenient to study the simplest gauge system, i.e. the 3-D Z 2 gauge model. It 
is dual to the 3-D Ising model where one could use non-local cluster algorithms for 
numerical simulations which yield a drastic reduction of the statistical error compared 
with the old Metropolis or Heat Bath methods. Direct measurements of expectation 
values of the Wilson loops with this new technique have not yet been performed, 
however there is a related observable that has been evaluated in this way. Let us see 
briefly the method. 

Denote by o\ = ±1 the link variables of the Z 2 gauge system and by a = ±1 the 
site variables of the dual Ising model. One can relate the expectation value of the 
plaquette P = OxG^OzO^ in terms of dual variables as follows 



with G\ = <Ti<t 2 , where / is the link of the dual lattice orthogonal to the plaquette P 
and punching the center of the plaquette; (3 is related to (3 by 



More generally, denoting by E an arbitrary surface with boundary <9E, we have 





(3 = — log (tanh (3) 



(9) 




(10) 



2.0 - 




l.O 1.5 2.0 2.5 3.0 3.5 

Fig. 1. Finite size effects in fluid interfaces. 

where the product is over all the dual links cutting the surface E. In the dual Ising 
model with periodic b.c. in x and y directions one can also realize a sort of Wilson 
loop which has no analogue in the gauge system, by choosing as E a toroidal surface 
wrapped around the periodic directions ( for instance a plane orthogonal to the z axis). 
Then (Ehes ex p( — 2/35";)} denotes the interface partition function. This observable 
is particularly well suited to study the shape dependent effects, because there are 
no complications due to the boundary sources. The size of this periodic Wilson 
loop is four times larger than the maximal ordinary loop one can measure in the 
corresponding gauge system. The functional form of Eq. 5 is now modified because 
the perimeter term is absent, of course, and the the gaussian term is slightly different 
with respect to Eq. 6 owing to the periodic boundary conditions. One gets 

q(x,y) = log (\/^smT\r](T)\' 2 } , r = i— . (11) 
v ' y 

A set of numerical simulations on the Ising model 14 and on the three-state Potts 
model 15 gave a nice, accurate confirmation of this behaviour for a wide range of r. 

In particular, the ratio R = Z(L\, L 2 )J Z(L } L) of interfacial partition functions 
for systems with the same area (LiL 2 = L 2 ) } but different shape, is only a function 
of L1/L2 and is directly expressible in terms of Eq. 11 with no free parameters. In 
Fig.l a comparison is made between such theoretical prediction and a set of numerical 
simulations in the Ising model at (3 = .2275 (black dots). 

3. Width of the Colour Flux Tube 

3.1. Thickness of the Flux Tube as an Effect of String Fluctuations 

One of the simplest and general consequences of the existence of a thin flux tube 
connecting a quark pair in the confining phase is that this flux tube fluctuates. The 



effect of these quantum fluctuations is to produce an effective square width of the 
flux tube which grows logarithmically with the interquark distance. 

Such a logarithmic behaviour has been predicted many years ago by Liischer, 
Miinster and Weisz 16 in the framework of the effective string picture of the gauge 
systems. In the next subsection we shall rederive and refine this universal law by 
directly using some exact results on two dimensional free field theory in a finite box. 

Let us begin by defining the observable we want to discuss. The square width (or 
the square gyration radius) of the flux tube generated by a planar Wilson loop W(C) 
is defined as the sum of the mean square deviation of the transverse coordinates of 
the underlying string, i. e. 

1 D ~ 2 r / x2 
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where T> is the planar domain bounded by C, A its area A = j v d 2 ^ and hf' M is the 
transverse coordinate of the center of mass of the flux tube, given by 



h ? M A 
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\j v d 2 ih l {( u ( 2 ) . (13) 



Combining Eq. 12 with Eq. 13 we get 

where Gi(£,£') is the Green function of the field hi defined as 

G i (te') = -(h i (Z)h i (Z')) . (15) 
a 

According to the discussion of §2, we expect that the 2D field theory describing 
the dynamics of the flux tube flows, for a domain T> big enough, to a massless gaussian 
limit, where the Green function fulfills the free field equation 

-AG(e, = ^-0 • (16) 

For simplicity we omitted the transverse index i. It is convenient to adopt complex 
coordinates by defining z = £i + ^2; then the solution G O0 (z } z') of the above equation 
in the infinite plane has the well known form 

Gooiz.z') = -77- log \z - z'\ + c , (17) 

Z7T 

where c is an arbitrary constant. Note that for z ~ z' the free Green function is 
ultraviolet divergent , then Eq. 14 is ill-defined and needs regularization. Using for 
instance the point-splitting method one gets at once for a square Wilson loop L X L 

w 2 = -L log (L/R c ) , (18) 
Zira 



where the UV cut-off has been absorbed in the definition of the scale R c . 

3.2. Green Functions in a Finite Box 

Actually the above calculation has a flaw, because the theory is only defined in a 
finite box, while it is used the propagator of the infinite plane. The Green function 
on a box has a very different form which depends on the shape of the box and on the 
choice of the boundary conditions. This fact might in principle modify the functional 
form of the width. Actually we shall see that this is not the case. 

The problem of finding the free Green function on an arbitrary, simply connected 
region T> of the plane with fixed boundary conditions has been already solved by the 
mathematicians of the last century. Accordingly, it can be reformulated in a way 
where z and z' play an apparently asymmetric role: find a real function f z i(z,T>) 
which is harmonic in the punctured set T> \ {z 1 } , vanishes at the boundary dT> and 
diverges logarithmically as — log \z — z'\ for z — > z 1 '. Such a problem can be solved 
once it is found a conformal mapping z — > I of T> onto the unit circle |£| = 1 which 
maps z' into the origin I = 0. Denoting by £ z >(z) the analytic function providing us 
with such a mapping, it is immediate to verify that 

= log |4'(*) | , (19) 

and that the function 

G v (z,z') = ^-f z ,(z,V) (20) 

Z7T 

satisfies Eq. 16. For, note that G-p is the real part of an analytic function, then 
AG-p = for z ^ z' ; since z' is mapped into the origin I = 0, then G-p diverges 
logarithmically for z — > z' and is normalized as required by the delta function; finally, 
being the boundary dT> mapped into the unit circumference, G-p vanishes at this 
boundary as it should. 

For our purposes, the relevant property of the conformal mapping I is that one can 
perform an arbitrary scale transformation z — > Xz without destroying the conformal 
character of I. More precisely we can write 

f z ,{z,V) =h z ,{\z,V x ) , (21) 

where T>\ denotes the scaled domain. A direct consequence is that it is always possible, 
for any finite T> } to fix the area of the scaled domain T>\ to an arbitrary value, say 
1, without changing the Green function. It follows that the integration of the finite 
part G(z } z') in Eq. 14 cannot depend on the size of the domain T> but only on its 
shape. Surprisingly enough, the logarithmic growth of the square width w 2 comes 
from the UV divergent part. Let us see why. According to the point spitting method, 
we regularize Cr(£, £) with the replacement 



(22) 



where e is the ultraviolet cut-off. Combining the logarithmic divergence of G with its 
scaling property expressed in Eq. 21, we can write 



G(i\e + £)~-^log(e/L) + c , (23) 

Z7T 

where L is a typical linear dimension of the domain T>. Inserting such an expression 
in Eq. 14 we get the logarithmic law of Eq. 18 where now R c is a calculable function 
of the shape of the domain. 

The physical meaning of such a behaviour is now clear: the dynamics of the flux 
tube is described by a truly free-field theory only at large distances; the cut-off e sets 
up the ultraviolet scale R c below which the free-field approximation breaks down. 
Obviously this scale cannot depend on the infrared scale L, then a variation of L 
cannot be balanced by a variation of e. In other terms the ultraviolet divergence 
breaks the scale invariance of the problem and then a variation of L gives rise to an 
observable effect. 

The volunterous reader may check the above arguments in a particularly simple 
example, given by the Green function on a disc of radius L 



Gdisc(z,z') = -77- log 
Z7T 
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(24) 



For a rectangular box of size L\ X L 2 , which is the most interesting case for the lattice 
gauge models, we have 



G LuL2 (z,z r ) = -T^-log 

Z7T 
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(z - z')a(z + z') 



a(z — z')a(z + z') 



(25) 



where (t(z) is the Weierstrass sigma function for the rectangle of sides 2L\ and 2L 2 . 
Inserting this expression in Eq. 14 it turns out that the scale R c does not depend very 
much on the modulus r = iL 2 / L\ for not too elongated rectangles. 

3.3. Numerical Simulations 

Though the logarithmic growth of the squared width of the flux tube is the most 
important and model-independent quantum effect predicted by the effective string 
description, it has not yet been observed until now. This is due to two main reasons. 
First, this effect is an infrared phenomenon that can be seen, as we have explained, 
only for very large Wilson loops, which are at the boundary of the sizes reached by 
the present numerical simulations. Secondly, a logarithmic growth can be checked 
numerically with reasonable confidence only if one can deal with a set of very precise 
data in a wide range of the loop sizes, which is again a condition which the current 
lattice simulations for SU (2) and SU (3) gauge theories have not yet attained. 

It has been recently observed 18 that in the simplest gauge system, i.e. the Z 2 3D 
gauge model, it is nowadays possible to overcome these difficulties by exploiting the 
duality relation between this model and the ordinary 3D Ising model explained in 
§2.2 . 



Fig. 2. Density of a 7j 2 flux tube evaluated at a distance of two lattice spacings from a Wilson loop 
24 X 24 sitting on the xy plane. Only the upper left quadrant is reported 



Indeed, using the one-to-one mapping of physical observables of the gauge sys- 
tem into the corresponding spin observables it is possible to replace the ordinary 
Metropolis or heath-bath method with a non-local cluster algorithm (for instance 
the Swendsen and Wang method). This allows us to probe the structure of the flux 
tube with an unprecedented accuracy. The procedure is the following. A Wilson 
loop W(C) is realized in the spin lattice by frustrating all the links cutting a given 
surface E bounded by C. These frustrated links modify the vacuum state so that 
the vacuum expectation value (P)w of the plaquette, or better its spin analogue, is 
no longer translational invariant, but becomes a function of its relative position with 
respect to W(L\, L 2 ) and is related to the expectation value in the ordinary vacuum 
by 

(P)w = (W(L u L 2 )P)/(W(L u L 2 )) . (26) 

The difference between the expectation value of the plaquette in the vacuum modified 
by the presence of W(L\, L 2 ) and the expectation value in the ordinary vacuum can 
be considered as a measure of the density of the flux tube. Choosing for instance as 
a probe a plaquette P|| parallel to the plane of the Wilson loop we can write 

p\\(x,y,z) = {P\\) w -{P) . (27) 

Other orientations of the plaquette give approximately the same distribution * The 
flux density p\\ can be viewed, up to an obvious normalization, as the probability for 

*This seems not to be the case in four dimensional systems 19,20 
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Fig. 3. Squared width of the flux tube in units of sigma. 

unit of volume of finding the fluctuating string at the point of coordinates (x } y } z). 
The main difference between the flux density in Z 2 gauge model and that in theories 
with a continuous gauge group is that in the discrete case there is no self-energy 
peak associated to the quark lines ; as a consequence, the formation of the flux tube 
can be seen more neatly. Starting from a point on the perimeter of the Wilson loop 
an moving inwards, one finds that the flux density increases rapidly in a few lattice 
spacings and then it reaches a stable plateau. The same feature can be found in any 
plane parallel to the Wilson loop, with a value of the plateau which is only a function 
of the distance from the plane of the loop as Fig. 2 shows. 

Then, assuming that the Wilson loop is located in the plane z = 0, the mean 
square width of the flux tube can be defined as 

w 2 = J dx 3 z 2 p\\(x,y,z) I J d x 3 p\\(x , y , z) , (28) 

where the integration volume is bounded by the plateau described above. We per- 
formed two sets of simulations. The hrst was at (3 = 0.7516 corresponding to 
era 2 = 0.0107(1) . We considered six different squared Wilson loops of sides ranging 
from L = 16 to L = 40 and fitted the data to the two-parameter formula 

w 2 = a log(L) + b . (29) 

We got b = — 17.4±4.0 and a = 14.4 zb 1 .2. The theoretical value a t h of the parameter 
a, fixed by Eq. 18 to be a t h = , matches nicely with this value, indeed we have 
a t h = 14.8 . The second set of simulations was made at (3 = 0.7460 corresponding to 
era 2 = 0.0189(1) with squared Wilson loops of sides ranging from L = 15 to L = 60. 
Fits to Eq. 29 gave b = —5.8 ± 2.0 and a = 7.7 ± 0.6 while the theoretical value is 
a th = 8.42 . 



^The high precision of this value is a consequence of the use of the non-local cluster algorithm. 



An internal consistency check of these numerical data comes from a comparison 
between Eq. 18 and Eq. 29. We get that the physical adimensional quantity \foR c is 
expressed in terms of a and b as 

^ = ^tL . (30) 

V Zira 

Using the fitted values of the parameters, we get \foR c = 0.35 ± 0.11 at (3 = 0.7516 
and \foR c = 0.31 ± 0.09 at (3 = 0.7560 in good agreement with scaling. 

One of the important features of Eq. 18 is that it can be written in a universal form 
by expressing all the dimensional quantities w, L and R c in units of \fo . Accordingly 
we report in Fig. 3 data from different /3's and and also recent very accurate data 21 
for the mean square width of fluid interfaces in the dual Ising model (black dots). 
The straight line represents the logarithmic fit to our data at (3 = 0.7460 (rhombs). 
Within the statistical accuracy these data clearly support a logarithmic widening of 
the flux tube in a range of quark separation L of more than two orders of magnitude, 
starting about L m i n \fo ~ 1.7. 

In order to compare these results with analogous data for other gauge groups, we 
can, with an abuse of language, express \fo in the same physical units of the QCD 
string . Then the crossover to the infrared logarithmic behaviour observed in Fig. 3 
corresponds to L m i n = .75 fm while the maximal probed elongation corresponds to 
more than 100 fm. 

The data on 4D SU (2) flux tube 19 ' 20 cover now a distance up to 2 fm, but are still 
affected by strong systematic errors and are compatible also with a constant width 
for separations larger than 1 fm. They are scattered mainly above the straight line 
of Fig. 3, suggesting a value of the UV scale a bit smaller than the value of the Z 2 
system, while the crossover distance is approximately the same. 

Such a similarity of the cross-over distance L m i n in so different gauge systems has 
a plausible explanation. According to Eq. 2, L m i n should be controlled by the first 
non gaussian term in the a' expansion the effective string action. This term does 
not depend on the nature of the gauge group, so it should produce universal effects. 
Actually it has been recently 22 shown that this term describes accurately some finite 
size effects in fluid interfaces which cannot explained in the free field limit. 
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